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Abstract. Let A be a fixed complex matrix and let u, v be two vectors. The 
eigenvalues of matrices A-\-ruv T (r E R) form a system of intersecting curves. 
The dependence of the intersections on the vectors u, v is studied. 

Introduction 

The motivation for this paper is the following numerical experiment. Take a 
matrix A e C™ xn and nonzero vectors u, v 6 C™ and plot the set 

(0.1) {a(A + ruv T ) :t eR}. 

It is well known that above set consists of a finite number of curves, that intersect 
only in a finite number of points. However, it appears that for w, v G C" chosen 
randomly from a continuous distribution on C ra there are no intersection points ex- 
cept, possibly, the of spectrum of A. Furthermore, for all r € R\ {0} all eigenvalues 
of A + ruv T , that are not eigenvalue of A, are simple. A typical case for A = Js(0) 
is shown of Figure [TJ note that the only intersection of the eigenvalue curves is at 
£ &{A). Since it appears that the intersection points outside <j(A) are multiple 
eigenvalues of A + tuv t (cf. Proposition 12. 2f ii)). we will be also interested in a 
problem of existence of multiple eigenvalues of A + tuv t for some r e C. 

Some light on the phenomenon of lack of double eigenvalues in the numerical 
simulations is put by the following marvelous result of Hormander and Melin [6]. 
Let the Jordan canonical form of the matrix A be 

r kj 

A — (BO-''' 

3 = 1 i=l 

where the Jordan blocks J nj , (Xj) corresponding to each eigenvalue Xj (j — 1, . . . r) 
are in decreasing order, i.e. n^i > nj.2 > • • • > fij^j ■ Then for generic u and v (i.e. 
for all u and v except a 'small' set, see Preliminaries) the Jordan form of A + uv T 
is the following 

r kj I 

A + uv T S J n}A (Xj) © JiHh, 

j = l i=2 h=l 

where fih ^ Hh 1 f° r h ^ h'. In other words, for each eigenvalue Xj (j = 1, . . . r) only 
the largest chain in the Jordan structure is destroyed and there appears a structure 
of simple eigenvalues instead. 



Key words and phrases. Perturbations, eigenvalues. 

The second author gratefully acknowledges the assistance of the Polish Minstery of Higher 
Education and Science grant NN201 546438. 

1 




FIGURE 1. Eigenvalues of B(t) = Js(O) + tuv t for t > in red, 
for r < in blue, see also Remark 14.21 

The behavior of eigenvalues of A + tuv t as functions of r for small values of r 
is also well known, see, e.g., [TU1IIH] and [U HI EH 031 E3 ■ Namely, for small values 
of | x and for generic u and v for each j = 1, . . . r there are rtj i simple eigenvalues 
fJ*j,k{T~), k = 1, ■ • • , n^i of A + tuv t in a punctured neighborhood of Xj, and they 
are given by 

(0.2) Mi , fc (r) = A, +r 1 ^ • (c,) 1 ^- 1 • exp (^) + 0(r 2 ^), 

where the number Cj can be expressed explicitly in terms of A,u and w; see |14) . 
Proposition 1. That is, the eigenvalues Hj t ki T ) are approximately given by the roots 
of the polynomial equation 

(0.3) (fi ~ Aj)"^' 1 — t ■ Cj, j \....r. 

However, neither the Hormander-Mellin result nor the above small r asymptotic of 
eigenvalues does not explain the lack of crossing of eigenvalue curves that appears 
in numerical simulations. The purpose of the present paper is to show that this 
behavior is indeed 'generic' although the notion of genericity will have some different 
shades. 

For historical reasons let us mention two works prior to the Hormander-Mellin 
paper, in [17] the invariant factors of a one-dimensional perturbation are considered 
and in [9] the perturbation theory for normal matrices is developed. The result 
by Hormander-Mellin lay dormant for about a decade before being rediscovered 
independently by Dopico and Moro [5] and Savchenko [TU [15]. Since that time 
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the interest in topic has grown up, see e.g. I12j for an alternative proof using 
ideas from systems theory and for perturbation theory for structured matrices. 
Although the results presented below concern a similar matter the reasonings are 
independent of the previous work and the content of the paper is self-contained. 
The main outcome are Theorems 13 . 1 1 |4~T1 15.11 16. II and 16.21 First four of them allow 
the parameter r to be complex, while in the last one we return to the real parameter 
r. This collection gives a complete description of the generic behavior of the set in 

1. PRELIMINARIES 

In this section, we gather some known results which will be the basis for our 
further investigation. An important technique used in this paper is the resultant. 
Let 

<Zi(A) = a ni A ni + ■ • • + a , <z 2 (A) = 6„ 2 A" 2 + • • ■ + b 

be two complex polynomials. By S(qi, q 2 ) we denote the Sylwester resultant matrix 
of qi and q 2 : 



(1.1) S(q u q 2 ) 















«() 



&o 









«o 




bo 



b„ 2 








Of) 






b 



f£>(rai+ri 2 ) X (ni+n 2 ) 



It is well known that q\ and q 2 have a common root if and only if det S(qi, q 2 
Let A G C nx ™ and let u,v G C™. Occasionally we will use the notation 

iv T , T G C, 



0. 



B(t) = vU 

remembering, nevertheless, that we are interested in the (it, independence of the 
spectral structure of B{t). Recall that an eigenvalue Ao of B G C nXTl is called 
non- derogatory if dimkcr(_B — Ao) = 1. The following result may be found in |14) . 
Lemma 5, for completeness sake we include a proof. 

Lemma 1.1. Let A G C nxn and let u,v G C n . Then for all t G C \ {0} all 

eigenvalues of B{t) that are not eigenvalues of A are non-derogatory. 

Proof. Let A G ct(B(t)) \ a {A) and let r ^ 0. Using the fact that rank (X + Y) < 
rank A + rankK for any compatible matrices X, Y we obtain 

n = rank (A — Ao) < rank (A + tuv t — Ao) + rank (tuv t ) — 

= rank(S(r) -A ) + l, 

which shows that rank(B(r) — Ao) > n — 1. Hence dimker(i?(T) — Ao) = 1 and so 
Ao is a non-derogatory eigenvalue of B(t). □ 

Following [11] we say that a subset il of C™ is generic if is not empty and the 
complement C" \ fl is contained in a (complex) algebraic set which is not C". In 
such case C" \ is nowhere dense and of 2n-dimensional Lebesgue measure zero. 
We use the phrase for generic v G C" as an abbreviation of: 'there exist a generic 
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f2 C C™ such that for all v E f2'. Our main results, except Theorem \6.2\ have the 
following form: 

Let A £ C" xn . Then for generic u and v 

which should be read formally as 

For every A € <C nxn there exists a generic subset fl ofC 2n , possibly depen- 
dent on A, such that for (u,v) € fL... 

Most of our reasoning are independent of a choice of basis. Let T be an invertible 
matrix. Then 

T(A + tuv t )T- 1 = TAT- 1 + t(Tm)(w t T" 1 ). 

In consequence, the Jordan structures of the matrices A + tuv t and TAT^ 1 + 
t(Tu)(v t T~ 1 ) are identical. In other words the transformation 

(1.2) (A, u, v T ) h-> (TAT-\Tu, v t T- 1 ) 

preserves the spectral structure of B{t) for all t£1. Let Ta be the transformation 
of A to its Jordan canonical form, that is 



(1.3) 



3 = 1 i=l 



where Jfe(A) denotes the Jordan block of size k with the diagonal entries equal A 
and the entries on the first upper -diagonal equal one and 



(1.4) 



n>j,i > n.j-,2 > • • • > n j:kj , j = 1, . . . , r. 



We will describe now a special instance of the transformation T that consists of 
two steps, i.e. T = T v Ta- Let Ta be as above, next we decompose u' = Tau and 
v' T = v T Tj according to the Jordan form of A' as follows: 



(1.5) 



and 



(1.6) 



We put 



*3,l 
4,2 



3,1 



l j,i,l 



"3,i,nj,t 



^ = 00Toep( Wj>i ), 
i=i «=i 

where by Toep (w) we denote the fc x k upper-triangular Toeplitz matrix whose first 
row is given by w G C k . Obviously T v commutes with A. Now note that for generic 
v one has 



(1.7) 



v 'i,i.i ¥= o 



1 , . . . , k j , j — 1 , . 
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which implies that T v is invertible, consequently T V A'T V 1 = A'. Furthermore, 



(1.8) 



= v T v has the following form 

"T 




l,...,kj, j 



The triplet (TAT' 1 , Tu, v T T- 1 ), where T = T V T A , will be called the Brunovsky 
form of (A, u,v T ), cf. [2J. Note the following simple lemma, that will allow us to 
reduce the problem of genericity in u and v to a problem of genericity in u with a 
fixed v. 

Lemma 1.2. I/Qq is a generic subset of C" then the set 

{(u,v) G C 2 ™ : T„ «s invertible ,T v Tau £ f^o} 
is a generic subset of C 2 ™ . 

2. The characteristic polynomial of S(t) 

The present section contains the basic tools used in the paper. Namely, we in- 
troduce the polynomial p uv and provide a formula for the characteristic polynomial 
ofS(r). 

The minimal polynomial of A will be denoted by m(A). Everywhere in the paper 
(|1.3|) and (|1.4p are silently assumed, consequently one has 

r 

m(X) = l[(\-X j r^. 

i=i 

p„„(A) = m(A) • v T (A - A)~ 1 u. 

Note that p u „ is invariant on the transformation (jl.2l) . Transforming A to it Jordan 
form we easily see that p uv is a polynomial of degree at most degm — 1. The 
following Lemma plays an essential role in the further reasoning. 

Lemma 2.1. For generic u and v the polynomial p uv is of degree degm — 1 and 
has no double roots and no common roots with m. 



(2.1) 

We also put 
(2.2) 



Proof. Using Lemma 11.21 and the fact that p uv is invariant on the transformation 
(|1.2j) we may assume that A is in the Brunovsky canonical form and treat v as 
fixed. For simplicity consider the case when A consists of one Jordan block only, 
i.e. 



Jn(Ai), v 





"1" 















V = 




, u = 














Thenm(A) = (A — Ai) 11 and 

(A -A)" 1 =Toep([(A-A 1 )- 1 ,(A- 
Consequently, 

Puv {\) =u 1 (X-X 1 ) n - 1 +■■■ 



Ai; 



■ , (A — Ai)~"] T ). 



u n -i(X - Ai) + u n 
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Hence, the generic assumption u\ ^ implies that degp U v = degm — 1. Further 
on, the generic assumption u„ ^ implies that p u and m do not have common 
roots. To prove that for generic u the polynomial p uv has simple roots only let us 
consider the Sylwester resultant matrix S(p uv ,p' uv ). Note that det S(p uv ,p' uv ) is a 
nonzero polynomial in u. Hence, the equation det S(j> uv ,p' uv ) = defines a proper 
algebraic subset of C™. 

The general case follows by similar arguments from the equation 

r kj 

p uv (X) = m(A) • J2J2 v Ja x - JnjA^T 1 ^^- 

j=l i=l 



a 



We put 



i=l j=2 y ' 



with the convention FJ* := 1- We also define the family of polynomials p UVtT by 
(2.3) Puv,tW = m(X) (A), rel. 

Proposition 2.2. Let A G C nxn , then the following statements hold. 

(i) For every u, v G C™, t G C the characteristic polynomial of A + tuv t equals 

Q ' PtLV,T • 

(ii) For every u, v G C", Ti,T2 G C m£/i n 7^ T2 owe ftas 

<t(A + TxUV T ) cr(A + t 2 uv T ) C 

(iii) for generic u and v and all t G C \ {0} t/iere are exactly deg m, counting 
algebraic multiplicities, eigenvalues of A + tuv t that are not eigenvalues of 
A. 

Point (iii) shows that the only crossings of the eigenvalue curves in (|0.1[) are the 
multiple eigenvalues of A + ruw T for some r G R. 

Proof, (i) For any u, v G C", r G C and A G C \ cr(A) we have (cf. |14j . Lemma 1) 
det(A - {A + tuv t )) = det ((A - A)(I - (A - A)" 1 ™/)) 
= det(A - A) det(I - (A - A)- 1 tuv t ) 
= dct(A - A)(l - tv t (X - A^u) 

Dividing both sides by q and emploing (|2.2[) we obtain 



/n ,n det A- (A + tuv t )) 

(2-4) i ^ = m(X)-T Puv (X), 

which finishes the proof of (i). 

(ii) Assume that Ao G a {A + r 1 uu T ) Pi a(A + t 2 uv t ) with t\ ^ r 2 . By (i) Ao is 
either a root of q, or a common root of the polynomials p uv ^ Tl and p U v,r 2 - I n the 
former case Ao clearly belongs to cr(A), in the latter case Ao is a root of (ri — T2)p uv 
and consequently of m. Hence, Ao G cr(^4) as well. 

(iii) By Lemma l2.ll for generic u and v and all r G C \ {0} the polynomials p uv . T 
and m do not have common roots and consequently q is the greatest common divisor 
of the characteristic polynomials of A and A + ruv T . Hence, for generic u and v 
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the roots of p uv , T are precisely the eigenvalues of B(t) which are not eigenvalues 
of A. Since the degp U v,r = degm, there are exactly degm,, counting algebraic 
multiplicities, eigenvalues of A + tuv t which are not eigenvalues of A. 

□ 

Note that by Lemma [1.11 for each r ^ the eigenvalues in ct(B(t) \ <j(A) are 
non-derogatory. However, the proposition above does not say, that for each t ^ 
the eigenvalues in <j(B(t) \ a(A) are simple. Obviously, for a fixed value of r and 
generic u and v the eigenvalues in o~(B(t) \ a (A) are simple, as follows from the 
Hormander-Mellin result, but this is a weaker statement. 

3. The Jordan structure of A + tuv t at the eigenvalues of A. 

The theorem below shows that the Jordan structure of B(t) at the eigenvalues 
of A is constant for all r ^ 0. The technique of the proof was used in pj] to reprove 
the Hormander-Mellin result . 

Theorem 3.1. Let A G C nx ™ and let ([O]), (LI} be the Jordan form of A. Then 
for generic u and v and all r G C \ {0} the sizes of the Jordan blocks of A + tuv t 
corresponding to the eigenvalue Xj are > ■ • ■ > n j.kj > f or J = 1 , . . - , r . 

Proof. Using the transformation (|1.2[) we can assume that A is in the Brunovsky 
canonical form. Denote by ejj (J — 1, . . . , r, I = 1, . . . , nj.i + nj.2 + • • • + nj,kj) 
the vector with one on the l-th position in the j-th block and zeros elsewhere. 
Then the following sequences are Jordan chains of A + tuv t corresponding to the 
eigenvalue Xj (j = 1, . . . r): 

e J,l ~ e j,rij,l+nj,2 + l> ■ ■ ■ i e j,rij,3 ~ e J,™3,l+"j,2+™3,3 i 

(3.1) . 

e jM ~ e 3)%',lH hrij.fcj -l + l i ■ • • J e j.n jikj ~ e i,%-,iH hnfc^ - 1 +rij,k j • 

Hence, we see that for generic u and v there are Jordan chains of A+tuv t of lengths 
n j.2 > • ■ • > n j.kj corresponding to the eigenvalue Xj. (Obviously, if kj = 1 then Xj 
is not an eigenvalue of A + tuv t ). By Proposition ^. 2l the dimension of the algebraic 
eigenspace corresponding to <j(B(t)) \ <j{A) is degm = n^i + • • • + n r ,i. Hence, 
none of the Jordan chains in (|3.ip can be extended and the proof is finished. □ 

4. The large t asymptotic of eigenvalues of B(t). 

In this section it is shown that the eigenvalues of B[t) that are not eigenvalues 
of A tend with r — > oo to the roots of the polynomial p uv , except one eigenvalue 
that goes to infinity. This behavior is again generic in u and v. 

Theorem 4.1. Let A G C" xn . Then for generic u, v G C n there exist differentiable 
functions 

: {t G C : \t\ > T } -> C, 
with I — deg m and some tq > 0, such that 

(i) <7(5(t)) \ o-(A) = { Mj -(r) : j = 1, ...,/} /or |r| > r ; 

(ii) //j ^ fi r for j, j 1 = 1, j ^ j'; 

(iii) yCti(r), . . . ,/i;_i(t) iend wit/i |t| — > oo £o i/ie ^ — 1 roots of the polynomial 

PltV ; 
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(iv) Hi(t)/t — > v T u with r — > oo. 

The theorem says, in other words, that as r goes to oo the eigenvalues of B(t) 
which are not eigenvalues of A are simple, exactly I — 1 of them are approximate the 
roots of p uv and one goes to infinity, asymptotically along the ray in the complex 
plane going from zero through the number v T u. 



Proof. By Lemma 12.1 1 there are I — 1 simple roots of the polynomil p uv , let us denote 
them by Ai, . . . , A;_i. Let e > be such that the closed discs 

C J (e) = {AeC: |A-A,| <e}, j = 1, . . .1 - 1 

do not intersect. Consider the polynomials 

<7r(A) = -m(A) - p Mt) (A), t > t 

T 

and observe that ^m(A) converges with |r| — > oo uniformly to zero on Uj=i Cj(e). 
By the Rouche theorem there is a tq > so that for |t| > tq the polynomial q T 
has exactly one simple root fJ,j(r) in each of the sets Cj(e), j = 1, . . . , I — 1. Hence, 
the root /i/(r) ^ Uj=i ^j( £ ) ^ s simple as well. By simplicity of the roots we get 
q' T (/ij(T)) ^ for j = 1,.../, |t| > To. Hence, by implicit function theorem the 
functions /xi(t), . . . , hi(t) are differentiable. Recalling that a{B(r)) \ a(A) consists 
by Proposition 12.21 precisely of the roots of (jv(A) finishes the proof of (i) and (ii). 
Letting e — > we obtain (iii). To prove (iv) note that 

V)WC±AWfi!M m, H >„. 



T I \T 



As t — > oo the matrix r 1 B{t) converges to the rank one matrix uv T and thus 
l JL i{ T )/ T converges to v T u. □ 

Remark 4.2. In Figure [1] the roots of the polynomial p uv are marked with black 
circles, and the asymptotic ray y — (v T u)x is the dashed line. 



5. Triple eigenvalues of B{t). 

In this section we show that for generic u, v there are no triple eigenvalues in 
ct(B(t)) \ a(A) for all r e C. In particular there are generically no triple crossings 
of the eigenvalue curves. 

Theorem 5.1. Let A e C" xn . Then for generic u,v £ C" and for all t € C the 
algebraic multiplicity of the eigenvalues of A + tuv t that are not eigenvalues of A 
is at most two. 

Proof. Suppose that u and v are such that for some r 6 C the matrix B(t) has an 
eigenvalue Ao ^ °~{A) of multiplicity at least three. Then by Lemma H. 1 1 B(t) has 
a Jordan block of size at least three at Ao- Consequently, by Proposition ^. 21 Ao is 
a triple root of p U v,n i-e. 

m(Ao) - Tp uv (\ Q ) = 0, 
m'(Ao) - rp^(Ao) = 0, 
m"(A o )-Tpl(A o ) = 0. 
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Solving for r from the first equation and substituting in the second and third we 
obtain 

m'(X )p uv (X ) - m(A )p' ut ,( A o) = 0, 
m"(\ Q )p uv (\ ) - m(\ )p" v (\ Q ) = 0. 

Let s be the greatest common divisor of m and m! . Since Ao does not belong to 
a (A), it is a common root of the polynomials 

Juv — Puv Puv> 
S S 

(5.2) g uv =m"p uv - mp\ 



uv ' 



Therefore, det S(f uv , g uv ) = 0. Summarizing, we showed so far that the set of all 
u and v for which there exists t G C such that the matrix B(t) has an eigenvalue 
Ao ^ v{A) of multiplicity at least three is contained in the set of all u,v £ C n such 
that det S(f uv , g uv ) — 0. Clearly det S(f uv , g uv ) is a polynomial in the coordinates 
of u and v. We show now, that it is a nonzero polynomial, i.e. that for some 
u,v the polynomials /„„, g uv do not have a common root, which will finish the 
proof. First consider the case degm = 1. Then for generic it, v the polynomial p uv 
is a constant nonzero polynomial and thus /„„ is a constant nonzero polynomial 
as well. Therefore, it does not have common roots with g uv . Now let us turn to 
the case degm > 1. Observe that for every b 6 C there exist such that 

Pu b v b W = A - b. Then 

fu bVb {\) = — (A)(A-6)--(A), 
s s 

SW(A) = m"(A)(A - b). 

Let fii,... ni-2 be the roots of m" . Note that —(//?) = implies —(fij) due 
to the definition of s. Therefore, one can find bo £ C \ a (A) such that 

TDj TTi TTlJ 

— in) ■ b o ¥= "(m?) - — (/^i) ■ Mji i = h-'-i - 2. 

Consequently, f Ub Vbg and g Ubo v bQ do not have a common root. □ 

Obviously, the result holds only generically. One can easily construct a matrix A 
and vectors u and v such that B(t) will have an eigenvalue of a given multiplicity 
for a given To. Namely, let Aq = Jfe(0) and let u,v be any two vectors for which 
A = Aq — tquv t has k different eigenvalues. Then A + tuv t — Jfe(0). 

6. Double eigenvalues of B(t) 

Theorem 6.1. Let A G C" xn . Then generic u, v G C™ there are at most 2 degm— 2 
values of the parameter r G C for which there exists an eigenvalue of A + tuv t of 
multiplicity at least two, which is not an eigenvalue of A. 

Proof. Note that for all r G K \ {0} the matrix B{t) has a double eigenvalue if and 
only if the polynomials p uv . T and p' uv T have a common zero, see Proposition 12.21 
Write m and p uv as 

l-i 2-i 
m(A) = A ; +^a J A- 7 , p ut ,(A) = ^ Pj A J . 
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Then the polynomials p uv , T and p' uv T (A) are given by 

l-l 

Puv,tW = A' + ^{a 3 -rp^X , 
3=0 

z-i 

P ' UViT (x) = ix 1 - 1 + J2j(*j -tp,)*' 1 - 

3=1 

Consider the Sylvester resultant matrix S(p UViT ,p' uv T ) G C' 2 ' _1 - |x ^ 2 '~ 1 - ) and let 

G(u,v,t) = det S(p uv , T ,p' UVT ). 

Then G(u,v,t) = if and only if there is an eigenvalue of -B(r) of multiplicity at 
least two, which is not an eigenvalue of A. Computing the determinant G(u,v,t) 
by development of (11.11) according to the first column (note that a ni = 1, b n2 = I), 
one sees that it is the sum of constant in r multiples of two determinants of size 
(21 — 2) x (21 — 2), the entries of which are linear polynomials in r, or constants. 
Using the fact that the determinant of a k x k matrix is a polynomial of degree 
k in the entries of the matrix, we see that G(u,v, ■) is a polynomial of degree at 
most 21 — 2 in the variable r. This means that for any A, u and v the polynomial 
G(u,v, •) has at most 21 — 2 zeros or is identically zero. However, by Theorem 14.11 
we already know that for generic u,v there exists r > such that for |r| > r 
the spectrum <t(B(t)) \ <r(A) consists of simple eigenvalues only and consequently 
G(u,v,t) ^ 0. Thus for generic u, v the polynomial G(u,v, ■) has at most 21 — 2 
roots and the theorem is proved. □ 

The last result of this paper considers the real parameter r. Together with 
Proposition I2.2f ii) it shows why the crossing of the eigenvalue curves in (jO.ip do 
not appear in numerical simulations, except possibly the crossings at <r(A). 

Theorem 6.2. Let A e C nxn and let V be the set of all pairs (u,v) G C 2 ™ for 
which there exists tGK such that A+tuv t has a double eigenvalue, which is not an 
eigenvalue of A. Then V is closed, with empty interior and has the An-dimensional 
Lebesgue measure zero. 

Proof. As in the proof of Theorem lG.ll we note that 

V = {(u,v) € C 2n : 3 remo} G(u,v,t) = 0} . 

Since the zeros of a polynomial depend continuously on its coefficients, the set is 
C \ V is open. To prove that V is of 4rt-dimensional Lebesque measure zero (and 
consequently has an empty interior) consider the set 

U := {(«,») € C 2 " : 3 AeC f uv (X) = £ (A) = 0} , 
where /„„ is defined as in (|5.1[) . Note that 

U = {(u, v) e C 2n : Bxec fuv(X) = guv(X) = 0} , 
where g uv is defined as in (|5.2p . Indeed, this follows from 

^ fuv — $9uv S f uv 

and from the fact that the polynomials s and /„„ do not have common roots. Hence, 
it follows from the proof of Theorem 16 . 1 1 that the set Uo is a proper algebraic subset 
ofC 2 ". 



EIGENVALUES OF RANK ONE PERTURBATIONS OF UNSTRUCTURED MATRICES 11 



Recall that by Lemma \2. II the set 

Ux = {(u,v) eC 2n :degp uv < I - l} , 

is also a proper algebraic subset of C". Observe that for (u,v) £ U\ one has 
deg/ u „ = k, where k := max{(r — — l),r(l — 2)}, / = degm and r = deg ^ is 
the number of eigenvalues of A. To see this let (u, v) £ U\. In the case (r— 1)(Z— 1) ^ 
r(l — 2) it is clear that degf uv = k. In the case when (r — 1)(7 — 1) = r(l — 2) 
note that although the degrees of both summands in (|5.1[) coincide, the leading 
coefficient does not cancel. Indeed, the leading coefficients of ^Puv and — p' uv are 
respectively la and (I — l)a, where a is the leading coefficient of p uv . 
Consequently, 

V := C 2 " \ (U U Ut) 
is an open and nonempty set. Note that for each (u, v) € Vo the function f uv 
has precisely k zeros Ai (u, v), . . . , Xk(u, v) and they are all not in cr(A). Since 
fuv(^j(u, v)) = and (u,v) ^ Uq, one has f' uv (\j(u,v)) ^ 0, j = 1, . . . , k. There- 
fore, by the implicit function theorem, the functions Ai(ti, v), . . . , A&(u, v) can be 
chosen as holomorphic functions on Vq. Note that 

k 

VCUoUUiUlJVj, 

3 = 1 

with 

Vj = {u 6 V : 3 reR \ {0} m(A 3 - («,«)) - T-p uu (A 3 -(u,v)) = 0} 
Pu«(A,-(M,t))) \ . 

Observe that the functions 

y 3 (ttj u) _ PMu,v)) = yT(x ^ y) A) _ lu e Cj ; , ; 

are holomorphic and nonconstant on every connected component of Vo- By the 
uniqueness principle each of the sets Vj (J = 1, . . . , k) is of 4n-dimensional Lebesgue 
measure zero. Hence, Vo, and in consequence V as well, are of 4n-dimensional 
Lebesgue measure zero. 

□ 

In the infinite dimensional case the function Q(z) = — ((A — A)~ 1 u, u) is a very 
useful tool for studying spectra of one dimensional perturbations of selfadjoint op- 
erators, or even more generally, spectra of finite dimensional selfadjoint extensions 
of symmetric operators. The key point is solving the equation Q{z) = — 1/r and 
as it can be seen this technique was a motivation for the proof above. This ap- 
proach can be found e.g. in [7] in the Hilbert space context and in [3j [4j [16] in the 
Pontryagin space setting. 

References 

[1] H. Baumgartel, Analytic Perturbation Theory for Matrices and Operators, Birkhauser, Basel 
1985. 

[2] P. Brunovsky, A classification of linear controllable systems, Kybernetika {Prague) 6 (1970), 
173-188. 

[3] V. Derkach, S. Hassi, and H.S.V. de Snoo, Operator models associated with Kac subclasses 
of generalized Nevanlinna functions, Methods of Functional Analysis and Topology 5 (1999), 
65-87. 



12 



A. CM. RAN AND M. WOJTYLAK 



[4] V.A. Derkach, S. Hassi, and H.S.V. de Snoo, Rank one perturbations in a Pontryagin space 

with one negative square, J. Fund. Anal. 188 (2002), 317-349. 
[5] F. Dopico and J. Moro, Low rank perturbation of Jordan structure, SIAM J. Matrix Anal. 

Appl. 25 (2003), 495-506. 
[6] L. Hormandcr, A. Melin, A remark on perturbations of compact operators, Math. Scand 75 

(1994), 255-262. 

[7] S. Hassi, A. Sandovici, H.S.V. de Snoo, and H. Winkler, One-dimensional perturbations, 
asymptotic expansions, and spectral gaps, Oper. Theory Adv. Appl. 188 (2008), 149-173. 

[8] T. Kato, Perturbation Theory for Linear Operators, Springer, New York 1966. 

[9] M. Krupnik, Changing the spectrum of an operator by perturbation, Lin. Alg. Appl. 167 
(1992), 113-118. 

[10] V.B. Lidskii, To perturbation theory of non— selfadjoint operators, Zh. Vychisl. Mat. i Mat. 

Fiz. (U.S.S.R. Comput. Math, and Math. Phys.) 6 (1966), 52-60. 
[11] C. Mehl, V. Mehrmann, A. CM. Ran, L. Rodman, Eigenvalue perturbation theory of classes of 

structured matrices under generic structured rank one perturbations, Lin. Alg. Appl. (2010), 

doi:10.1016/j.laa.2010.07.025. 
[12] C. Mehl, V. Mehrmann, A. CM. Ran and L. Rodman, Perturbation theory of selfadjoint 

matrices 2 and sign characteristics under generic structured rank one perturbations, Lin. 

Alg. Appl. (2010), doi:10.1016/j.laa.2010.04.008. 
[13] J. Moro, J.V. Burke, M.L. Overton, On the Lidskii- Vishik-Lyusternik perturbation theory 

for eigenvalues of matrices with arbitrary Jordan structure, SIAM J. Matrix Anal. Appl. 18 

(1997), 793-817. 

[14] S.V. Savchenko, On a generic change in the spectral properties under perturbation by an 
operator of rank one, [Russian] Mat. Zametki 74 (2003), 590-602; [English] Math. Notes 74 
(2003), 557-568. 

[15] S.V. Savchenko, On the change in the spectral properties of a matrix under a perturbation of 

a sufficiently low rank, [Russian] Funkcional. Anal, i Prilozhen. 38 (2004), 85-88; [English] 

Funct. Anal. Appl. 38 (2004), 69-71. 
[16] H.S.V. de Snoo, H. Winkler, M. Wojtylak, Zeros and poles of nonpositive type of Nevanlinna 

functions with one negative square, arXiv:1011.2081 
[17] R.C Thompson, Invariant factors under rank one perturbations, Canad. J. Math. 32 (1980), 

240-245. 

[18] M.I. Vishik, L.A. Lyusternik, Solutions of some perturbation problems in the case of matrices 
and self-adjoint and non-self-adjoint differential equations, Uspekhi Mat. Nauk (Russian 
Math. Surveys) 15 (1960), 3-80. 

[19] J.H. Wilkinson, The Algebraic Eigenvalue Problem, Oxford Univ. Press, Oxford, 1965. 

Afdeling Wiskunde, Faculteit der Exacte Wetenschappen, Vrije Universiteit Ams- 
terdam, De Boelelaan 1081a, 1081 HV Amsterdam, The Netherlands 

Instytut Matematyki, Wydzial Matematyki i Informatyki, Uniwersytet Jagiellonski, 
Krakow, ul. Lojasiewicza 6, 30-348 Krakow, Poland 



